Astronomy & Astrophysics manuscript no. ki09may astroph 


© ESO 2009 


November 3, 2009 





Letter to the Editor 



o 
o 

> 

o 
m 



O 

u 

O 



> 

o 



o 



X 



Exploring the inner region of Type 1 AGNs witli thie Keck 

interferometer 

Makoto Kishimoto^, Sebastian F. Honig^, Robert Antonucci^, Takayuki Kotani^, Richard Barvainis^, Konrad 

R. W. Tristram\ and Gerd Weigelt^ 

'Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany; e-mail: mkOmpifr-bonn.mpg . de. ^Physics 
Department, University of California, Santa Barbara, 93106, USA. ^ISAS, JAXA, 3-1-1 Yoshinodai, Sagamihara, Kanagawa, 229- 
8510 Japan. ''National Science Foundation, 4301 Wilson Boulevard, Arlington, VA 22230, USA. 



Submitted 20 October 2009 / Accepted 3 November 2009 



ABSTRACT 



The exploration of extragalactic objects with long-baseline interferometers in the near-infrared has been very limited. Here we re- 
port successful observations with the Keck interferometer at K-band (2.2 yum) for four Type 1 AGNs, namely NGC4151, Mrk231, 
NGC4051, and the QSO IRAS 13349+2438 at z=0.108. For the latter three objects, these are the first long-baseline interferometric 
measurements in the infrared. We detect high visibilities (V^ ~ 0.8-0.9) for all the four objects, including NGC4151 for which we 
confirm the high level measured by Swain et al. (2003). We marginally detect a decrease of with increasing baseline lengths 
for NGC4151, although over a very limited range, where the decrease and absolute are well fitted with a ring model of radius 
0.45±0.04 mas (0.039+0.003 pc). Strikingly, this matches independent radius measurements from optical-infrared reverberations that 
are thought to be probing the dust sublimation radius. We also show that the effective radius of the other objects, obtained from the 
same ring model, is either roughly equal to or slightly larger than the reverberation radius as a function of AGN luminosity. This 
suggests that we are indeed partially resolving the dust sublimation region. The ratio of the effective ring radius to the reverberation 
radius might also give us an approximate probe for the radial structure of the inner accreting material in each object. This should be 
scrutinized with further observations. 
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1. Introduction 

The exploration of extragalactic objects, or in particular. Active 
Galactic Nuclei (AGNs), with long-baseline interferometers in 
the near-infrared (near-IR) has been very limited. While the 
brightest Type 1 AGN NGC4151 and Ty pe 2 AGN NGC1068 
have been observed by Swain et al. (2003) and Witt kowski etal] 
( 12004) . respectively, further exploration has been hampered 
mainly by technical difficulties. Here we report successful ob- 
servations of four Type 1 AGNs with the Keck interferometer 
(KI) in the near-IR (K-band 2.2 fim). Type 1 AGNs are thought 
to give us a direct view of the innermost region of the putative 
dust torus as well as the central accretion disk, where the inter- 
esting effect of the latter should also be evaluated carefully. 



2. Keck interferometry 

2.1. Observations and data reduction 

We observed four AGNs listed in Table [T]and associated calibra- 
tors w ith the Keck Interferometer (Kl HColavita & WizinowichI 
l2003h on 2009 May 15 (UT). These four targets were chosen 
based on their bright optical magnitudes measured from the pre- 
imaging data obtained in April 2009 at Tiki Observatory (French 
Polynesia) and Silver Spring Observatory (USA) by N. Teamo, 
J. C. Pelle, and K. Levin. 

The KI combines the two beams from the two Keck 10 m 
telescopes which are separated by 85 m along the direction 
38° east of north. Adaptive Optics coiTection was implemented 



at each telescope, locking on the nucleus in the visible wave- 
lengths. The data were obtained with a fringe tracker rate of 200 
Hz operated at K-band, while the angle-tracking was performed 
at H-band. The data were first reduced with Kvis' to produce raw 
squared visibility (V^) data averaged over blocks of 5 sec each. 

Then those blocks with a phase RMS jitter larger than 0.8 ra- 
dian were excluded. For a given visit of each object, the blocks 
with an AO wavefront sensor flux smaller than the median by 
20% or more were also excluded. The rejected blocks were ~ 
6% of all the blocks, and they were generally outliers in the visi- 
bility measurements for a given visit, and quite often were of low 
fractional fringe-lock time. Then the wide-band side of the data 
were further reduced using wbCalib^ with the correction for the 
flux ratio between two telescope beams and the coiTection for a 
flux bias (slight systematic decrease of the KI's system visibility 
for lower injecte d flux^). The jit ter correction was applied with a 
coeflicient 0.04 (IColavitallT999l) . Then the blocks were averaged 
into scans over a few minutes each, with its error estimated as a 
standard deviation within a scan. 

FiglT] shows observed visibilities of all the targets and cal- 
ibrators (after the coiTections above), plotted against observing 
time. All the six calibrators used are expected to be unresolved 
by the KI at K-band (V^ > 0.999). Overafl, the system visi- 
bility (V^ys), as measured by these calibrators, was quite sta- 
ble over the night. The calibrators span over a relatively wide 
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range of brightness (see the legend in Fig[T), and one of them 
(HDl 11422) had approximately the same injected flux counts 
as those of NGC4151 and Mrk231. Based on the corrected vis- 
ibilities of these calibrators shown in Fig[T] the flux bias seems 
to have been taken out quite well, although there might still be 
some systematics left. The visibilities of the three faint calibra- 
tors (K > 9.1) tend to be slightly smaller than those for the other 
brighter ones, with the diff'erence of the means of the former and 
latter being ~0.9% of the means. Therefore we assign 0.01 as a 
possible systematic uncertainty in system visibility estimations. 

Note that for the faintest target NGC4051, the flux bias cor- 
rection was effectively used with a slight extrapolation (by ~1 
magnitude in K-band^). This should be checked with future 
fainter calibrator observations. 

For each target measurement, V^^ was estimated from these 
calibrator observations with wbCalib using its time and sky 
proximity weighting scheme, yielding V^y^ indicated as gray cir- 
cles in Fig[T] The final calibrated visibilities are shown in Fig|2l 
with the sampled uv points shown in the inset of Fig|2] 

2.2. Results 

Figl3] shows the final calibrated visibilities for NGC4151 as 
a function of projected baseline length, enlarged in the left 
panel of the insets . We confirm the visibility level observed by 
ISwain et"al ] (l2003h . The covered range of projected baselines is 
very limited (note also that the shortest possible for NGC4151 
is ~70 m due to the KI's delay line restriction). We see, how- 
ever, a marginal decrease of visibility over the increasing base- 
line. With the Spearman's rank correlation coefficient analysis, 
the confidence level is 98.4%, or 2.4 cr. The decrease and ab- 
solute level of visibilities are well fitted with a simple thin ring 
model (i.e. inner radius equal to outer) of radius 0.45 ± 0.04 mas 
(0.039 ± 0.003 pc; the error accounts for the systematic uncer- 
tainty in V^y^; Sec l2.1l ). If we convert each of the visibility mea- 
surements into a ring radius and plot it as a function of the PA 
of the projected baseline, we obtain the right panel of the insets 
in Figl3] Over the PA range covered, from ~10° to ~50°, we do 
not seem to see PA dependence of the radius. For all the other 
targets, the calibrated visibilities are shown in Table [T] together 
with baseline information and deduced ring radii. 



3. UKIRT imaging and Kl data corrections 

In order to obtain quasi-simultaneous flux measurements for the 
nuclear point-source in the four Type 1 AGNs (which are vari- 
able), their images were obtained with WFCAM on UKIRT in 
five broad-band filters (Table|2]i on 2009 June 17 (UT) under the 
UKIRT service programme. The seeing was ~1 arcsec. The pipe- 
line-reduced data were obtained through the WFCAM Science 
Archive. The wide field of view of the WFCAM gave simultane- 
ous measurements of PSF stars in each AGN field, while a 2x2 
micro stepping gave a good image sampling with an effective 
pixel size of 0.2 arcsec. We implemented two-dimensional (2D) 
fits for each image to accurately separate the PSF component 
from the underlyi ng host galaxy , follow ing the same procedure 
as described by K ishimoto et alj (l2007h . 

Table |2] lists the measured flux of the nuclear PSF compo- 
nent. The uncertainty of our nuclear PSF flux measurements is 
estimated as ~5%, based on the residual fluxes after the fits and 
the flux calibration uncertainty. In the K-band images of the two 
brighter objects, namely NGC4151 and Mrk231, the central sev- 
eral pixels seemed affected by non-linearity or saturation, so we 
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Fig. 1. Observed plotted against observing time. Gray dots 
are individual measurements for blocks of 5 sec each. Gray cir- 
cles are the estimated system at the time of target observa- 
tions. 
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Fig. 2. Calibrated plotted against observing time. The inset 
shows the sampled uv points for each target (north to the top, 
east to the left). 



implemented the 2D fits by masking the central ~0.4 arcsec ra- 
dius region. We estimate that the nuclear PSF flux is recovered 
within the same uncertainty of 5%, based on the results from the 
same masked fits on the other unsaturated images. Using the re- 
sults of the PSF-host decomposition, we also estimated the host 
galaxy flux fraction within the field-of-view of the KI which is 
~50 mas at K-band (FWHM; set by a single-mode fiber for the 
fringe tracker). The results are stated in Table |2] The obtained 
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Table 1. Properties of our targets and summary of the results of our KI observations on 15 May 2009 (UT). 



corr. 



scale 

(pc mas"') 



Eb-v'^ 
(mag) 



B 



P 
(m) 



PA 

(°) 



(mas) 



(pc) 



(pc) 



Av 
(mag) 



NGC4151 0.00414 0.086 0.028 85.0 12.3 0.862±0.018 0.45±0.04 0.039±0.003 

NGC4051 0.00309 0.064 0.013 79.4 41.9 0.861±0.026 0.51±0.05 0.032±0.003 

MRK231 0.0427 0.84 0.010 74.4 44.3 0.923±0.028 0.38±0.07 0.32±0.06 

IRAS13349+2438 0.109 2.0 0.012 85.0 34.5 0.869±0.016 0.44±0.03 0.88±0.05 



~0.044±0.011 
-0.011±0.()04 



1.3" 
0.93' 



" CMB corrected value from NED. * //q = 70 km s"' Mpc"', £1^ = 0.3, and = 0.7. Galactic reddening from lSchlegel et all <1998h . Projected baseline 
lengths. For objects with multiple measurements, data shown are at the longest baseline length, except for NGC405 1 where the one with a smaller error is 
shown. ^ Thin-ring radius (best-fit values if multiple data are available). ^ Mean and standard deviation of the reverberation measurements over 2001 to 2006 from 
iKoshida et aiU20^ for NGC4151. and over 2001 to 2003 from lSuganuma et all HmS ) for NGC4051. * iLacv et al J fT982h . ' IWills et all fT99l . 
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Fig. 3. Calibrated for NGC4151 as a function of projected 
baselines, enlarged in the left inset. The dotted line shows the 
best-fit visibility curve with a ring model of radius 0.45 mas. 
In the right inset, ring radii are plotted along the PA of each 
projected baseline. Note that the correction for the accretion disk 
and host galaxy contributions is not incorporated in this figure, 
and it does not change the ring radius significantly (see Table|2]i. 



small values show that the host galaxy contribution is only a 
very small part of the observed visibility departure from unity. 

FigH] shows the resulting spectral energy distribution (SED) 
of the PSF component in each target, after the correction for 
Galactic reddening. We also corrected for the reddening in the 
host galaxy for the objects which show large B aimer decrements 
in broad emission lines (Table |2l). Assuming that the PSF flux 
originates from the hot dust thermal emission nearly at the sub- 
limation temperature and from the central accretion disk (AD; 
thought to be directly seen in Type 1 inclinations), we esti- 
mate the flux fraction at K-band from the latter AD compo- 
nent. Here we fit the SED with a power-law spectrum of the 
form fy oc v^^^^ for the AD, plus a spectrum of a black-body 
form for the dust (the best-fit temperature was ~1300-1500 K). 
The AD flux fraction at K-band is estimated t o be as small as 
~0.2 (Table|2l), in agreement with the results bv lKishimoto et all 
( 120071) . This suggests that the high visibilities observed are not 
due to t he unresolved AD, as opposed to the preferred interpreta- 
tion bv lSwain et all ( |2003|) for NGC4151. The assumed near-IR 
AD spectral form is based on the recent study of near-IR po- 
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Fig. 4. Flux of the nuclear PSF component in WFCAM images 
derived from 2D fits. Fitted SEDs are shown in dotted lines (see 
text). 



larized flux spectra (iKishimoto et al.L l2008h . but also on various 
studies on AD spectral shapes in the optical/UV (summarized 
in Fig. 2 of Kishimoto et al. 2008). By assigning the uncertainty 
in the AD near-IR spectral index as 0.3, we also estimated the 
uncertainty of the K-band AD flux fraction (Table|2|. 

Finally, we corrected the observed visibilities for the host 
galaxy and the AD contributions, where the latter is assumed 
to remain unresolved by the KI. The corrected as well as 
corresponding thin ring radii are listed in Table |2l 



4. Interpretations and discussions 

We interpret here the high visibilities observed for all the four 
objects as an indication of partially resolving the inner bright- 
ness distribution of dust thermal emission. 

As we discussed in Secl3] the AD flux fraction at K-band 
is estimated to be small, as long as the assumed power-law AD 
spectrum in the near-IR, smoothly continuing from the optical, 
is at least roughly correct. In this case, the K-band emission is 
dominated by the dust emission, and it is reasonable to convert 
the observed visibility to a thin ring radius to obtain an approxi- 
mate effective radius of the dust brightness distribution for each 
object. (We have corrected for the unresolved AD contri- 
bution, but the correction is quite small; Table |2]) The derived 
ring radii are plotted in parsec in Fig|5] against UV luminos- 
ity L, here defined as a scaled V-band luminosity of 6 vfy{V) 
(Kishimoto et al., 2007). The V-band flux is extrapolated from 
the fitted flux at 0.8 fim (FiglH assuming an AD spectral shape 
of /y oc v" (based on spectral index studies referred to above). 
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Table 2. Point-source flux from UKIRT imaging data and the KI results corrected for host galaxy and unresolved AD component. 



name 






flux (mJy) 








magnitude 




host" 


AD 


V'- corrected 


i?,ing corr. 


fit" 




Z 


Y 


J 


H 


K 


J 


H 


K 


(%) 


fraction'' 




(pc) 


(pc) 


NGC4151 


39.7 


42.4 


59.4 


96.3 


173. 


11.0 


10.1 


8.9 


0.3±0.1 


0.17±0.06 


0.841±0.024 


0.041 ±0.004 


0.033 


NGC4051 


8.88 


9.98 


14.9 


24.8 


46.1 


12.6 


11.6 


10.4 


1.5±1.0 


0.14±0.05 


0.870±0.038 


0.032±0.005 


0.011 


Mrk231 


21.0 


22.0 


33.5 


71.6 


152. 


11.7 


10.4 


9.1 


0.5±0.1 


0.15±0.05 


0.921±0.033 


0.33±0.07 


0.33 


IRAS 13349 


7.75 


8.71 


12.6 


26.3 


62.3 


12.7 


11.5 


10.0 


0.4±0.1 


0.14±0.04 


0.856±0.021 


0.92±0.06 


0.49 



" Host galaxy flux fraction at K-band estimated for the KI's 50 mas FOV in the AO-corrected images. * AD flux fraction of the point source at K-band. ' Rj^ from 
UV luminosity using the fit bv lSuganuma et all 420061) . It has an uncertainty of a factor of ~ 1.5 based on the scatter of the fit (see text). 
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Fig. 5. Corrected ring radius derived for each KI target (squares), 
plotted against UV luminosity, or a scaled V-band luminosity 
(extrapolated from the WFCAM Z-band flux; see text). Also 
shown in gray plus signs are the rev erberation radii ag ainst the 
same scaled V-band luminosity (Sug anuma et al.L [20061 and ref- 
erences therein) and their fit (dotted line). 



We can directly compare these ring radii with another type 
of independent radius measurements Rj^, namely the light trav- 
eling distance for the time lag of the K-band flux variation from 
the UV/optical variation (Suganuma et al., 2006). These rever- 
beration radii are also plotted against the same scaled V-band 
luminosity in Fi gJS] They are known to be approximately pro- 
portional to L'^^ dSuganuma et al.ll2006l : the dotted line in Fig|5] 
shows their fit), and are likely to be probing the dust sublimation 
radius in each object. We first see that R^i^^ is roughly compara- 
ble to /?TK for all the four objects, and thus ^ring is roughly scaling 
also with l}^'^ . This approximate match suggests that the KI data 
are indeed partially resolving the dust sublimation region. 

With a closer look at Fig|5] we see that R^^^ is either roughly 
equal to or slightly larger than Rj^ (i.e. Rrmg/RrK ^ 1, up to a fac- 
tor of a few), though we have only four objects. This could be un- 
derstood if Rj^ is tracing a radius close to the innermost bound- 
ary radius of the dust distribution. It is known that the cross- 
correlation lag tends to trace an inner radius of the responding 
particles' distribution when the lag is determined from the peak 
in the cross-correlation function (e.g. Koratkar & Gaskell, 1991, 
and references therein), as is the case for the data used in the fit 
by Suganuma et al. On the other hand, Turing is an effective, aver- 
age radius over the radial dust brightness distribution in the K- 
band. When the radial distribution is very steep and compact, the 
ratio Riing/Rr^ would become very close to unity (such as seen 
in NGC4151 and Mrk231), while for a flatter, more extended 
distribution, Rnng/RrK would show a larger departure from unity. 



If our interpretation above is correct, the KI data would con- 
versely support the dust sublimation radius as probed by the re- 
verberation measurements that is smaller by a factor of about 
three than that inferred for typical ISM-size Graphite grains 
(BarvainiilHlTl 0.05 fim radius) for a given L (iKishimoto et al.L 
[2007 ). The small sublimation radius could be due to the possible 
dominance of large grains in the innermost region, since they can 
sustain at a much closer distance to the illuminating source for a 
given sublimation temperature. Altern atively, it could be due to 
an anisotropy of the AD radiation (see IKishimoto et aLll2007l for 
more details). 

In Fig(5| the reverberation and ring radii are shown essen- 
tially as a function of an instantan eous L at the time o f each 
corresponding radius measurement. iKoshida et all (l2009b . how- 
ever, recently have shown that is not exactly scaling with the 
instantaneous L varying in a given object. It might be that R^^ 
tends to give a dust sublimation radius that corresponds to a rel- 
atively long-term average of L. On the other hand, 7?^^ does show 
the L'^^ proportionality over a sample of objects. Thus, when we 
compare R^^ and Rnns, unless simultaneous measurements exist, 
we would have to allow for the uncertainty in TJ^k, as a function 
of instantaneous L, being the scatter in the L'^^ fit (~0.17 dex). 

If the AD spectrum does not have a power-law shape but 
rather has some red turn-over in the near-IR, though the near- 
IR p olarized flux spectrum argues against this (IKishimoto et al.L 
'2008), the AD flux fraction at K-band would become higher 
than we estimated here. However, even if the flux fraction is as 
large as 0.5, the corrected ring radius would become larger than 
stated in Table |2] by a factor of only ~1.3, resulting in no qual- 
itative change in our discussion. Future near-IR interferometry 
with much longer baselines can conclusively confirm that the 
visibility is decreasing as we inferred from the present KI data. 
We plan to advance our exploration with further interferometric 
measurements in the infrared. 
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